High-temperature expansion for steady nonequilibrium states in driven lattice gases 
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We develop a controlled high-temperature expansion for nonequilibrium steady states of the driven 
lattice gas. We represent the steady state as P{rj) oc e^W) - *^ and evaluate the lowest order 
contribution to the nonequilibrium effective interaction ^(n)- We see that, in dimensions d > 2, 
all models with nonsingular transition rates have the same ^(jy), which consists of a three-body 
effective interaction decaying like l/r 2d+1 . The models with the Metropolis rule show exceptional 
behavior. 
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To construct statistical mechanics that apply to 
nonequilibrium states has been among the most challeng- 
ing open problems in theoretical physics. So far it is not 
known how the desired probability measures for nonequi- 
librium states look like, or even whether the measures can 
be written in universal and compact forms. 

In a standard approach to nonequilibrium systems, one 
studies a concrete many-body problem with a given mi- 
croscopic dynamics, hoping to get some hints for possible 
universal statistical descriptions. Among the most well 
studied models is the driven lattice gas Ll [J . in which 
particles on a lattice are driven into one direction by an 
external electric field. For us, a major issue concern- 
ing that model is to elucidate universal aspects of the 
nonequilibrium steady state. 

Here, we develop a high-temperature expansion for the 
nonequilibrium steady state of the driven lattice gas. We 
represent the steady state as in ||S} with the nonequilib- 
rium effective interaction ^(n). We explicitly compute 
^(r]) in the lowest order of the expansion, without mak- 
ing any unjustified assumptions (such as the truncation 
of higher order correlations) . Although there were some 
attempts to develop expansions in the driven lattice gas 
0, 0, 01 , our approach to evaluate the effective interac- 
tion explicitly seems to be novel. Recently field theoretic 
perturbation method was developed for one-dimensional 
soft core lattice gas models driven by boundary condi- 
tions [||. 

We will find that, in the lowest order (0(/3 2 ) for d > 2, 
and 0(f3 2 ) or 0(/3 3 ) for d = 1), the effective interaction 
can have only two-body and three-body interactions as in 
Interestingly, these effective interactions generally 
are non-local and exhibit power law decay even though 
the Hamiltonian is completely short ranged. Such ob- 
servations may be of fundamental importance when one 
seeks for a universal description of nonequilibrium steady 
states. 

When the transition rate function <j>{h) is non-singular 
and the dimension is higher than one, our result of the 
expansion is rather satisfactory, suggesting a universal 
behavior of the nonequilibrium steady states. We see 
that in the lowest order, the nonequilibrium effective in- 



teraction ^(rf) consists only of an effective three-body 

(3) 

interaction tpx,y,z and is completely independent of the 
transition rule and the density. The effective interaction 

decays as in Ijl7(l and is summable, i.e., Yl y z IV^j/^l < 00 
in the infinite volume limit. 

The models with the Metropolis update rule, however, 
show quantitatively different behavior. This suggests 
that the choice of transition rule can be a serious problem 
in driven nonequilibrium systems, as was pointed out in 
0,0- Note that the Metropolis rule has been used widely 
in numerical computations of the driven lattice gases. 

Model: Let A C Z d be the d-dimensional L x • • • x L 
hypercubic lattice with periodic boundary conditions. 
We denote the sites as x = (xi, . . . , xj) £ A. A configu- 
ration of the system is described by r\ = (r/aOzeA, which 
is a collection of occupation variables r) x . We set n x — 1 
if x is occupied by a particle and r] x = if x is empty. 
Throughout the present paper, we fix the lattice size L 
and the particle number N — X^eA 7 ^- The interaction 
between the particles is described by the Ising Hamilto- 
nian H(rj) = -(J/tyJ^ix^VxVy When we write J2( x ,y)' 
we sum over all pairs of sites with \x — y\ = 1, double- 
counting x, y and y, x. 

The stochastic dynamics of the driven lattice gas is 
determined by specifying the transition rates. For neigh- 
boring x, y, we set (as usual) 

c(x ^y; V ) = &$lv<l>(E{rf«)-H(Ti) + E(xx-yi)), (1) 
which is the rate for the particle to hope from x to y in the 
configuration r\. Here, 6xX y — T] x (l — rj y ) is 1 only when 
the hop is possible, rf v is the configuration obtained by 
switching r\ x and r\ y in the original configuration 77, and E 
is the electric field in the x% direction. (We have absorbed 
the inverse temperature [3 into J and E.) We use the 
periodic boundary conditions so that x\ — y\ — 0, ±1. 
The function 4>(h) satisfies <fi(h) = e~ h (f)(—h) so that the 
local detailed balance condition 



x;v x ' v ) (2) 



c( x -> y-n) = ^(■n)-BW")-E(x 1 - Vl ) c f :y 

holds. We also require 4>(h) = 1 + 0(h). Then the above 
condition implies 4>'(0) = — 1/2 provided that <fi(h) is dif- 
ferentiate at h = 0. The standard choices are i) the 



2 



exponential rule with <p(h) = e~' l//2 , ii) the heat bath 
(or Kawasaki) rule with <p(h) = 2/(1 4- e' 1 ), and iii) the 
Metropolis rule with <p(h) — 1 for h < and = e - ' 1 
for ft > 0. 

Steady state: The steady state distribution P(r/) is 
the unique solution of 



-P(v) <x ^y; V ) + P{if y ) c{y -> i; = 0, 



or, equivalently, 



(3) 



p(rf,v) 



P(n) 



for any 7/ (with the fixed N = j] x ) . Let us write the 
steady state distribution as 



P(rj) 



Z(J,E) 



(5) 



where ^(n) is the "effective interaction" which represent 
a nonequilibrium correction to the measure. Then Q 
becomes 

= ]T c(y - z; rf ) ^W-^") - c(z - y; 77), (6) 

which is the starting point of our perturbation theory. 

Since we have ^(rf) — when E = or J = 0, we 
shall develop an expansion in and J, and (here) find 
the lowest order contribution to the effective interaction 
^(77). (We can as well think that the inverse temperature 
appears in the model as J = (3J and E = f3E, and one 
expands in (3.) We first note that the left-hand side of 
© is written as {1 + 0(E) + O(J)} A$(rj), where 

A*(r?) = ]T 6it v {*(v x ' v ) - *(»?)}. (7) 

Let us denote the right-hand side of © as —Q(i~i). 
From (0, we have 

Q{V) = E C 1 - e^" 1- * ) cfc -> y; 77). (8) 

We first assume that <j)(h) is differentiable at h — 0. By 
substituting (JTJ into JSJ and expanding in E and J to 
the second order, we get Q(rj) ~ Qgcn(v) with 

Qgen(ry) = §£ 41, (yi - *i) - (9) 

(a:>l/) 

where the result is independent of the specific choice 
of 4>(h). For the Metropolis rule, which has non- 
differentiable (f>(h), the result is different. We get Q{rf) ~ 



<3mp(t7) with 

Q M p(r?) = -£^ {H( V ) - H(rf*)}, (10) 

fail/) 

xi>yi 

provided that E > (2d— 1)| J\. The formula is much more 
complicated for \E\ < 2(d - l)\J\. 

We have thus found that the lowest order contribution 
to the effective interaction ^(77) is determined by the 
Poisson-like equation 



A*(t?) = -Q(tt), 



(11) 



where we set Q(rj) = Q ge n(v) or Q(v) = Qmp(t?) depend- 
ing on the transition rules. 

Decomposition: The equation (|llf) looks intractable 
since it is formulated in the huge space of all possible 
configurations 77. Fortunately, there are beautiful de- 
composition properties which enable us to determine the 
solution rather easily. 

We start from the Laplacian J7J). For Ac A, let rj A — 
UxeA Vx- We expand *(t?) as *(t?) = J2aca ^ a V A with 
real coefficients ipA- Note that an arbitrary function of 
77 can be expanded in this way. From the translation 
invariance and the particle number conservation, we can 
assume i)ja = when \A\ = 1. (We denote by \A\ the 
number of sites in A.) Now note that (rf' v ) A — j] A:c y 
for any 77 and A C A, where A Xjy is obtained from A 
by switching the occupation status at x and y. (If, for 
example, x € A and y ^ A then A Xiy = A\{x} U {y}-) 
Then from the linearity of the Laplacian J7J , we see that 

Av f fo) = EE^ ^ (^•* - v A ) 

(x,y) A 

= EE(K,-«^i" 

= ^(AVm)^, (12) 



where 

Aip A = El {V^Ux}^} ~ ^A>- (13) 

y^A 

We then treat the lowest order charges 10 and 

CIJJ. Note that 5 x %{H(r)) - H(rj x ' y )} = rj x (l - 

VyjJxYlz (\z-x\=i,zjty) T lz ~ (\z-y\=i,z^x) Vz}- Since 
the charge Q(rf) is a linear combination of these terms, 
it can be expanded as 

Q(V) = J2 $1 + IxizVxVyVz, (14) 

x,y£A x,y,z£A 

(2) (3) 

with suitable coefficients q x , y and q x ,y,z- Note that when 

(2) (3) ' ' 

evaluating the "charges" q x ,y, Qx,y,z, we only need to con- 
sider configurations rj with two or three particles. This 
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means that the calculations are (in principle) elementary. 
In what follows we only present the results of charge cal- 
culations, leaving details to @- 

By substituting these decompositions into the Poisson 
equation Ijlll) for ^(n), we find that ipA = unless \A\ = 
2 or 3, and tpA with \A\ = 2, 3 are determined by the 
(tractable) Poisson equations 

av41, = ( 15 ) 

Here Laplacian A is defined by ((T^|l by identifying {x, y} 
or {x,y,z\ as the subset A. Equivalently they may 
be regarded as the standard Laplacian Q on the lat- 
tices A 2 = {{x,y)\x,y £ A, x / y} C Z 2d or A 3 = 
{(a;, y, z) | x,y, z G A, x ^ y, y ^ z, z ^ x} C Z 3d . Two 
sites in A 2 (resp. A3) are regarded to be neighboring 
when their Euclidean distance in Z 2d (resp. Z 3d ) is equal 
to 1. 

The effective interaction ^(n) is exactly written as 



^IzVxVyVz, (16) 



x,y ,z£ A 



in the lowest order. It is remarkable that the effective 
interaction (at least in this order) is independent of the 
particle number. We stress that we are not expanding in 
the density. We shall now determine ip x , y and ipx,y,z by 
solving the Poisson equations (|15fl . 

General models in d > 2: We first concentrate on 

general models with (j>(h) being differentiable at ft, = 0. 

(2) 

As for the two-body charge, we find that q x . y = for 

(2) 

all a;, y. This means ip X} y = 0. We have no two-body 
correction to the measure in the lowest order. 

(3) 

The three-body charge q x , y , z , on the other hand, is 
nonvanishing. Let e\ = (1, 0, . . . , 0), and W be the set of 
the 2(d — 1) unit vectors of Z d other than ±e\. One has 

q x 3 l± ei x+ s — ±JE for any x £ A and any 5 £ W . Other 

(31 

nonvanishing Qx,y,z are determined by requiring it to be 
symmetric in x, y, and z. The rest are vanishing. To get a 

(3) 

feeling for the behavior of the solution ip x ,y,z of the Pois- 
son equation (JTSJ, we make an orthogonal transformation 
and introduce a new coordinate by X = (2x — y — z)/ 
Y = (y — z)/V2, and Z = (x + y + z)/\/3. Apparently 

the charge distribution depends only on X and Y, and so 

(3) (3) (3) 

does the solution ip x ,y,z- By denoting ip x ,y,z = fxy ' we 



get Atp%L 



(3) 

Qx yi wn i cn is the Poisson equation in the 

(3) 

2c?-dimensional space. Here the charge distribution q X Y 
is still complicated, but when projected onto the (Xi, Yi) 
plane (where X\ and Y\ are the first coordinates of X 
and Y , respectively), we see that the plus charges are lo- 
cated at (a/273,0), (-l/V6,±l/\/2) and minus charges 
at (-a/273, 0), (1/V6,±1/a/2) . Close to the origin of 
the 2d dimensional space, there is a hexa-pole parallel to 
the (Xl,Yi) plane. Roughly speaking the corresponding 

solution tp^y behaves like {l/ r 2d - 2 }'" - l/r 2d+1 . This 
estimate can be made into a reliable asymptotic estimate 



|8| by going back to the original 3d-dimensional lattice, 
and treating properly the "mirror charge" induced at the 
defect sites. The result is 



y.z 



cJE 



XJxVxWx 



V 



\w\ 



2 \d+2 ■ 



(17) 



where c is a constant depending only on the dimension, 
and we have set U — 2x — y — z, V — 2y — z — x, and W = 
2z — x — y. We note that quantity in the denominator can 
be written as |f7| 2 + |^| 2 + |M/| 2 = 2{\x-y\ 2 +\y-z\ 2 +\z- 
x\ 2 }. It is notable that the \/ r 2d+1 decay of l(T7|l implies 

2Z y z \ipx,y,z\ < const, with an L independent constant. 
We recall that such a summable (effective) interaction 
should be regarded as "healthy" in the light of the well 
established theory of the infinite volume Gibbs state ^(| • 

To summarize, the lowest order contribution to the ef- 
fective interaction ^(n) is independent of the transition 
rule provided that d > 2 and 4>'(0) exists. Moreover ^(n) 
in this order consists only of the summable three-body 
interaction with the asymptotic behavior l|17l) . The ro- 
bustness of the result may suggest that we have captured 
a universal aspect of the nonequilibrium steady state. 

It should be noted, however, that the higher order con- 
tributions to ^(n) may be more complicated and ma" 
become rule-dependent. In particular the results in 



may 
in 

( 2 ) 



suggest that the two-body effective interaction is 
nonvanishing and exhibits a l/r d power law decay as in 
(|18fl in the order of E 2 J for the heat-bath rule, while 
there is no l/r d decay for the exponential rule. 

Metropolis rule in d > 2: In models with the 
Metropolis rule (with E > (2d — 1)| J|), the situation 
becomes completely different. By the translation invari- 
ance, we have q x 2 y = q x 2 l yl where gi 2 ' = EJ if z = ±ei, 

q[ 2) = -EJ if z = ±2ei or z = ±ei + S, q {2) = 2EJ if 

(2) 

z = S, for any u £ U , and q z = otherwise. Again, 

by the translation invariance the corresponding effective 

(2) (2) (2) 
interaction can be written as ip x ,y = <Px—y> w here tp z 



is determined by Aip 1 ' 2 ' 1 = —qt } ■ (Here A is the Lapla- 
cian on A\{(0, . . . , 0)}.) Since q^ corresponds to the 
charge distribution of a quadrupole in <i-dimension, the 
resulting field tp^ is expected to decay like l/|z| d for 
large \z\. Again, this can be made into a precise estimate 
[3, and we have 



^( 2 ) 



tpW ^c'EJ\z\-^{(d-l) Zl 7 



i=2 



(18) 



for any d > 2 where c' is a constants which depend only 
on the dimension. This power law decaying effective in- 
teraction can be regarded as the origin of the l/r d long 
range two-point correlation 3> S S 0. 01 known to exist 
in the driven lattice gas and related models. Note that 
the interaction decays so slowly that the sum X^eA \ v^,v I 
is divergent as L — > 00, which implies (and indeed proves 
[TTI| ) that the infinite volume limit of the steady state 
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is not Gibbsian. Thus, in a stark contrast with general 
models with nonsingular <j>(h), the steady state measure 
for Metropolis dynamics has a significant nonequilibrium 
correction in the two-body sector already in the lowest 
order. 

Curiously enough, in the lowest order of the perturba- 
tion, the behavior of the three-body effective interaction 
Vx,y,z is exactly the same as the general models that we 
have discussed. We still do not know if this is a mere 
accidental coincidence or an indication of a deeper uni- 
versality in the structure of nonequilibrium steady states. 

Models in d = 1: The models in d = 1 show excep- 
tional behavior, which requires extra estimates 0. 

Suppose that 4>{h) is twice differentiable at h = 0. In 
this case nonvanishing contribution to ^(r/) is found only 
in the third order. The formula 10 for the charge is thus 
no longer useful. By going back to the original definition 
©, we get the following results. In the two-body ef- 

(2) 

fective interaction ipx,y> there is a short-range correction 
proportional to {20" (0) - {\/2)}E 2 J. The three-body 

effective interaction ip%*,y,z exhibits a l/r 3 long range cor- 
relation as in l|17fl . and its magnitude is proportional to 
(j)"(0)EJ 2 . Unlike in the models with d > 2, the low- 
est order correction is already very much rule dependent. 
Note in particular that <£"(0) = in the heath bath rule. 

In the Appendix 4 of [lj , one-dimensional exactly solv- 
able models of driven lattice gases are introduced, where 
the nonequilibrium steady states coincide with the cor- 
responding equilibrium states. In the language of the 
present work, these models turn out to posses vanishing 
charge ©. We therefore recover the known results that 
the models have no nonequilibrium corrections. 

As for the one-dimensional model with the Metropolis 
rule, there is nonvanishing contribution to ^(rj) in the 

(2) 

order EJ . We find a short range correction to ip x ,y, and 
no contribution to v4 3 y,z. 

Discussions: We have developed a high-temperature 
expansion for the driven lattice gas, and evaluated the 



effective interaction ^(rj) to the lowest order. Unfor- 
tunately, we still do not know how to evaluate higher 
order contributions in a systematic manner (although a 
brute force calculation seems always possible). We ex- 
pect to get higher body effective interactions as we go 
to the higher orders in the perturbation. We note that 
the expansion is much easier in the models with soft-core 
interactions, where we can take into account the effect of 
the field E without expanding it |8(|. It is also possible 
[U to write down equations for *(r?) in lattice gas models 
driven by chemical potential difference at the boundaries, 
where we expect the l/r d ~ 2 as in 0] to be universal [||. 

We have observed, as in 0,0, that the nonequilibrium 
steady states sometime show qualitative update-rule de- 
pendence. It is especially sever in one-dimension, and in 
the Metropolis rule in any dimensions. Such rule depen- 
dence can be a serious problem especially when one is in- 
terested in understanding universal statistical properties 
of nonequilibrium steady states. It is certainly desirable 
to find a criterion that determines the rule which yields 
physically meaningful results. 

We found that ^(r?) in d > 2 is independent of the 
rule provided that 4>{h) is differentiable. This strongly 
suggests that the results thus obtained are universal and 
physically meaningful ^3 ■ As we go to the higher orders 
in the expansion, however, even the models with differen- 
tiable 4>(h) are expected to show strong rule dependence 
Q • We still do not know of any way to determine the 
"right" transition rule. It is our conjecture 00 that the 
models with the exponential transition rule are free from 
the 1 /r d long range correlation, and have chance of pos- 
sessing Gibbsian infinite volume steady states. Although 
there seems to be no physical reasons to prefer models 
with Gibbsian steady state, it would be quite interesting 
if such a criterion plays a role in the future development 
of nonequilibrium statistical mechanics. 
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